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Abstract. The relativistic hydrodynamic system of equations for a perfect fluid 
obeying a causal equation of state is hyperbolic [T]. In this report, we derive the 
conditions for this system to be convex in terms of the fundamental derivative of the 
equation of state [5] • The classical limit is recovered. 
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Recent applications in several fields of relativistic astrophysics demand the use of 

realistic equations of state (EOS) beyond the polytropic and ideal gas EOS. In this work, 

we report on the characterization of convexity of general relativistic hydrodynamics 

(GRHD) and, more specifically, on the conditions imposed by thermodynamics. 
Ou 3 f [va 

Let — — I — = be a scalar non-linear conservation law (in one spatial 

at ox 

dimension, ID) for the conserved variable u : R x R — y R, u — u(x,t), f : R — > 
^) / = fi u ) being a non-linear function of u (flux). The behaviour of f(u) has 
important consequences on the behaviour of the solution u(x,t) of the conservation 
law itself [3J. A crucial property is the monotonicity of the characteristic speed \{u) = 
df/du = f'(u). A convex (concave) flux satisfies A' (it) = dX/du = f"(u) > (< 0), 
Viz. If X(u) has extrema, i.e., 3u such that X'(u) = dX/du = f"(u) = 0, then the 
conservation law describes a non- convex, non- concave flux. 

Let us now consider a non-linear, hyperbolic system of conservation laws (HSCL) 
of p equations in ID 

^ + ^ = (1) 

dt + 3x ' 1 ' 

where u : R x R — > MP, u = u(x, t) is the vector of conserved variables, and 
/ : MP — y MP, f = f(u), a non-linear function of u, the vector of fluxes. Let A Q (u) be 
the (real) eigenvalues and r Q (u) (a = l,2,...,p) the corresponding right eigenvectors of 

the Jacobian matrix — — . 

ou 

A characteristic field C a of system satisfying 

C a : ^ = A Q (2) 

is said to be genuinely nonlinear or linearly degenerate if, respectively, 

V a := V u X a -r a ^0, (3) 

V a := V u A a ■ r a = 0, (4) 

for all u, where V u A a is the gradient of A a (u) in the space of conserved variables and 
the dot stands for the inner product in MP. 

Conditions d3J) and (j4j), introduced by Lax [1] are the extension to systems of the 
above property on monotonicity of the characteristic velocity described for the scalar 
case [3J, and the HSCL is said to be convex (comprising both convex and concave fluxes 
in the scalar case) if all its characteristic fields are either genuinely non-linear or linearly 
degenerate. Additionally, in a non-convex system, non-convexity is associated with those 
states u for which one V a is zero and changes sign in a neighbourhood of u. 

In classical fluid dynamics, the convexity of the system is determined by the EOS 
[2l [3J and, more specifically, by the so-called fundamental derivative, Q (see its definition 
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and properties in, e.g., [2]) 



Q: 




dp 
dV s 



(5) 



p being the pressure, V := 1/p the specific volume (p is the rest-mass density) and s 
the specific entropy. An alternative expression for Q is [2] 



where c SjC ; := y dp/dp\ s is the classical (i.e., non-relativistic) definition of the speed of 
sound. 

It is important to note, however, that in fluid dynamics, one speaks of convexity 
when the thermodynamics forces rarefaction waves to be expansivdjj. In practice, this 
corresponds to impose the sign of Q and, as we shall see below and in the Appendix, 
also the one of the V a associated to the genuinely nonlinear fields. The fundamental 
derivative measures the convexity of the isentropes in the p — V plane and if Q > then 
the isentropes in the p — V plane are convex and the rarefaction waves are expansive. 
Then, one speaks of convex EOS and the result is a convex system, with characteristic 
fields which are either genuinely nonlinear or linearly degenerate, and definite signs for 
the V a associated with the genuinely nonlinear fields. Conversely, a non-convex EOS 
(in which Q < or has not a definite sign) leads to a non-convex flow dynamics^. 

The main goal of this report is to characterize the convexity of GRHD by examining 
the products V a defined in ([3]) and (Tj0) (Lax's criterion). From this analysis, the 
convexity will be described in terms of a new condition on Q that generalizes the 
classical Menikoff-Plohr's result [2]. The Appendix establishes the same condition on 
Q but obtained following Menikoff-Plohr's approach applied to self-similar relativistic 
flows in ID, and discusses the connection of this result with the choice of the signs of 
V a associated with the genuinely nonlinear fields. Our results generalize those outlined 
in ®. 

The evolution of a relativistic fluid is governed by a set of conservation laws, namely 
the conservation of rest mass, V ■ J = 0, and the conservation of energy-momentum, 
V ■ T = (V- stands for the four- divergence). For a perfect fluid, the components 
of the rest-mass current, J, and the energy-momentum tensor, T, are J M = pu 11 , 
and T^ v = phu^u u + pg^, respectively, h being the specific enthalpy, defined by 
h = 1 + e+p/ p, where e is the specific internal energy, is the four- velocity of the fluid 
and g^y defines the metric of the spacetime M. where the fluid evolves. Greek (Latin) 

% A rarefaction wave is said to be expansive when fluid elements decrease their pressure and density 
as they go through it. 

§ However, note that for more general HSCLs, there can be non-convex states even for a convex EOS, 
as in the case of classical ideal magnetohydrodynamics [5]. 



G = 1 + 



d In c Sid 
<9 In p 



s 



(6) 
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indices run from to 3 (1 to 3), or, alternatively, they stand for general coordinates 
{t, x, y, z) ({x, y, z}). The geometrized system of units (c — G—1) and the summation 
convention over repeated indices are used. 

An EOS p = p(p, e) closes the system. Accordingly, the (relativistic) sound 

speed c s := J dp/de\ s , e := p(l + e), satisfies he 2 = x + with x '■= dp/d p\ £ and 

k := dp/de\ p . 

In [7], the equations of GRHD were written as a hyperbolic system of conservation 
laws within the framework of the {3 + 1} formalism (see, e.g., [8]). According to this 
formalism, the metric is split into the objects a {lapse), (3 l (shift) and 7^, keeping the 
line element in the form: 

ds 2 = -(a 2 - ^)dt 2 + 2(3 i dx i dt + 7ii dxW. (7) 

If n is a unit timelike vector field normal to the spacelike hypersurfaces T, t (t = 
const.), then, by definition of a and /3\ d t = an + (3 l di, with n • di = 0, Vi. Observers, 
Oe, at rest in the slice i.e., those having n as four-velocity [Eulerian observers), 
measure the following velocity of the fluid: 

where := —(u^n^) = an 1 , the Lorentz factor, satisfies = (1 — v 2 )^ 1 / 2 , with 
v 2 = fjf* = 7i/u J ). 

The set of conserved variables gathers those quantities which are directly measured 
by Oe, i-e., the rest-mass density, D, the momentum density in the j-direction, Sj, and 
the total energy density, E. In terms of the primitive variables w = (p, Uj, e), they are 



D = pW, Sj = phW 2 Vj, E = phW 2 - p. (9) 

With the above definitions, the fundamental first-order, flux-conservative system ruling 
the evolution of flows in a given spacetime (GRHD equations) reads 



1 /<9 v / 7 F°(w) <9^/=#F i (w) 



yj—g V 9x° dx % 

where 



s(w), (10) 



u := F°(w) = (A^r) (11) 

is the vector of conserved variables, where r := E — D (total energy density excluding 
the rest-mass one), 

F- ( w) = ( D („• - g , s , („• - g + ^, x („• - g + ^) ( i 2) 

are the fluxes in each spatial direction, and the corresponding sources s(w) are 
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s w 



(o. r u) . ° - ) • (13) 



In the previous expressions, g := det(g^ v ) and 7 := det(7jj), and are such that 

The three 5x5- Jacobian matrices B l associated with system ( ITUj) are 
d¥ i 

The full spectral decomposition {A n ,r n }W (n = 1,...,5) of the above matrices B % can 
be found in [9j [J|. In terms of the primitive variables w, we can define the following 
matrices 

A> = „*™, A°= 9 -^. (15) 
aw aw 

It has been shown in [101 that: 



a) The matrices A i and & satisfy & = A\A )- 1 

b) If {A*, r*}w (n = l ; 5) are, respectively, the eigenvalues and eigenvectors of the 
characteristic eigenvalue problem corresponding to the system written in quasi- 
linear form, (A 1 — X*A°)r* = 0, then the following relations are satisfied 

{A n , r n }« = {A; , A <}« (n = 1, 5). (16) 

On the other hand, 

Proposition 1. Given an arbitrary 3-vector Q, the general eigenvalue problem 
(QA' — A*4°)r* = has the following eigenvalues and right eigenvectors: 

( ■ ft\ 

• Aq = a£j [v J \, degenerate (triple) eigenvalue associated with the material 



waves; 



a 



associated with the acoustic waves; 

• r*o,m = {"^iCmiiX)} where m (= 1,2,3) identifies the three independent right 
eigenvectors associated to X , and Cmi (^Vn = 1,2,3^ stands for the zero vector and 
the two unitary 3-vector s forming an orthonormal basis with Q/ WCjUi 

• r* ± = (p«A±,«(Cy -A±)((i-A ± v l ),-p X A ± -phW 2 (( j vi -A±) 2 ) ; with 
A± = -W 2 (A± - Cy) 2 - (QC j ) + AL and A± = (A± + Qp)/a. 

■ 

The right-eigenvectors published in [7] correspond to those spacetimes having a diagonal metric. 
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The spectral decomposition in the j spatial direction is recovered by taking 

The special relativistic limit (Minkowski 



d = o~l and, therefore, C = Y 8 J k 
spacetime) of the above expressions is recovered (in Cartesian coordinates) just doing 



a 



1,^ = 0,7^ 



;i7) 



Finally, taking into account that 

V u A Q (u) • r Q (u) ^ <=^> V w A a (w) • r* Q (w) ^ 0, 
we can formulate the following proposition: 

Proposition 2. The quantity V±, in GRHD, can be given in terms of three factors, one 
coming from the particular renormalization procedure followed to get the eigenvectors, 
TZ±, a kinematical one, K,±, and a purely thermodynamical one, T : 

v ± = n ± K ± r, (is) 

where 1Z± = apnA±. The kinematical term is: 

, W 2 



5 



with 



5 = W[{\ 
The thermodynamical term is: 



(I - c 2 s )(Qv j ) 2 ] 



(19) 



21 1/2 



T(p, e = — + — — + - 1 - cj 

op p l oe p 



(20) 



The term 1Z±K,± has the same sign for all the states (u) of the system. It only tends 
to zero when a — > 0, W — > oo (ultrarelativistic flow speeds), or c s — > 1 (ultrarelativistic 
thermodynamics) . 

The thermodynamical term in V±, T, can be written, in a more compact way, as a 
function of the fundamental derivative, Q (see Eq. (jSJ)) 

3 



(21) 



This formula has been derived considering: i) the definition of Q given in terms of the 
classical speed of sound (see Eq. (JH}), ii) the relationship between both definitions of 
the local speed of sound, c s c i 2 = he 2 , from which it is easy to show 



dine 



s,cl 



<9 hip 



,d 

<91nc= 



<91np 



(22) 



and iii) the thermodynamical partial derivative: 



dQ 


_ dQ 


p dQ 


dp 




e p 2 de 
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where Q stands for any thermo dynamical function of p, e. Hence, as far as the 

thermodynamical term is concerned, Eq. ( 121]) on one hand, and Eqs. ( 133]) . (134)) of the 

Appendix on the other, allows us to establish the following corollary: 

Corollary 1. Unlike classical fluid dynamics, which is convex for Q > 0, GRHD is 

3 

convex when the inequality Q > — c 2 is satisfied. 

Alternatively, in terms of the (relativistic) fundamental derivative Q, GRHD is 
convex when 

d 2 p 



~ 1 h 22 d(h/p) 2 

y : = ~o _ l 1 - c s) 



2 p v SJ dp 



d(h/p) 



> 0. (23) 



It can be easily seen that for causal (i.e., c 2 < 1) EOS, the derivative dp/d{h/ p)\ 



s' 



appearing in the denominator of the definition of G, must be negative. Hence, for 
causal EOS, the convexity condition expressed in (123"]) translates to d 2 p/d(h/p) 2 \ s > 0. 
Israel [H] proved that d 2 (h/ p) / dp 2 \ g > is a necessary condition for the stability 
of compressive shocks and the non-occurrence of rarefactive shocks. In the same 
context, Lichnerowicz [12] showed that the compressibility conditions d(h/ p) /dp\ s < 0, 
d(h/p)/ds\ p >0 and d 2 (h/ p) / dp 2 \ s > are sufficient conditions for a shock to be 
compressive. Lately, Thorne [13] (see also [11]) proved the same result as Lichnerowicz 
eliminating the second condition and only for weak shocks. Taking into account the 
nature of shocks as limiting simple waves we can conclude that our result, expressed in 
Corollary 1, comprises these early results on the characterization of compressive shocks. 

It would be interesting to probe the present result in numerical experiments and 
astrophysical applications using complex EOS. 
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Appendix. Characterizing the RHD convexity through Menikoff-Plohr's 
approach and the connection with Lax's criterion 

The SRHD (special relativistic hydrodynamics) equations for mass and momentum 
conservation in ID are respectively 

d( P W) d( P Wv) _ 

dt + dx _U ' 1 j 
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d{phW 2 v) d(phW 2 v 2 + p) 

di + die " ° (25) 
(where v stands for the flow velocity along the x direction, and the rest of quantities 
were already defined). 

Self-similar solutions of the above equations for isentropic flow (dp = he 2 , dp) are 
obtained by imposing that all variables are functions of £ = x/t. The following relations 
between the differentials are easily obtained: 

(v - £)dp + phc 2 [l + vW 2 (v - 0] dv = 0, (26) 
[1 + vW 2 (v - 0] dp + phW 2 (v - £)dv = 0. (27) 
Eliminating the differentials, the following condition has to be fulfilled 

(v - 0W 2 = Tc s [I + vW 2 {v - 0] , (28) 

which can be used to obtain the relation between the differentials dp and dv, 

dp = ±phW 2 c s dv. (29) 



Equation (1281 allows us to write £ as 

= (30) 
1 ± vc s 

Now, differentiating the previous expression and using the thermodynamic relation 
dh = dp/p, valid for isentropic flow, it is straightforward to show that 

1 — v 2 ~ dp , . 

(l±vc s y phc s 



where 

/ ) 1 > i / • 

" cl (32) 



Sine, 

1 + 



<91np 

which can be easily written in terms of the fundamental derivative, Q, using Eqs. ([6]) 
and (122) 

<? = S-^. (33) 

At this point, it is important to note that can be written as 

d 2 p 



lh 2 , 2 d(h/p) 2 



2p y s ' dp 



(34) 



d(h/p) 

in complete correspondence with the definition of Q, (JS}. 

Convexity (i.e., the fact that rarefaction fans are expansive) implies that the 
derivative d£ + /dp must be positive, and the derivative d£_/dp must be negative [2] 
f „ rci n get „ b e P o siti ve,o r6 >5,, , fc ^«,,, W <» md 
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d£_/dp > 0, rarefaction waves are compressive. If the sign of these derivatives is not 
defined and d^±/dp is positive for some states and negative for others the acoustic waves 
lose the genuinely non-linear character. 

All the above analysis can be done in a compact way by using the spectral 
decomposition of the ID SRHD system. Particularizing the spectral decomposition 
presented in Proposition 1 to this case, the corresponding eigenvalue problem in terms 
of the primitive variables w = (p, v, e) leads to the following eigenvalues 

A = v, A± = ^— ° s (35) 
1 ± vc s 

and eigenvectors 

r$ = (-«,0,x), (36) 

r ± = ( 772 IT! „2!,„2 J ' ( 37 ) 



he 2 ' phcgW 2 ' p 2 hc 2 s 

where the normalization of the eigenvectors has been chosen so that its first 
component, r* ± 1; fulfills the equation 

-¥- = dp (38) 
r ±,i 

for isentropic flows. 

Now, taking into account that, in a self-similar flow, the self-similar variable is 
£± = A± and the following relation between the variation of the variables across the flow 
and the components of the right eigenvectors holds, 

dp dv de 



±,1 '±,2 '±,3 



(39) 



the previous election allows us to write the differential of the eigenvalues A± as 

d\± = V w A±(w) • dw = V W A± • r*± dp. (40) 
Therefore, the sign of V W A± • determines the character of the flow. Convexity 

— * 

(expansive rarefaction fans) is obtained for a positive value of V W A+ ■ and a negative 
value of V W A_ ■ r^. 
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